**FULL TITLE** 

ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** 
**NAMES OF EDITORS** 



2 types of spicules "observed" in 3D realistic models 



Juan Martmez-Sykora 

Institute of Theoretical Astrophysics, University of Oslo, Norway 
Lockheed Martin Solar & Astrophysics Lab, Palo Alto, USA 

Abstract. 

Realistic numerical 3D models of the outer solar atmosphere show two 
different kind of spicule-like phenomena, as also observed on the solar limb. The 
numerical models are calculated using the Oslo Staggered Code (OSC) to solve 
the full MHD equations with non-grey and NLTE radiative transfer and thermal 
conduction along the magnetic field lines. The two types of spicules arise as a 
natural result of the dynamical evolution in the models. We discuss the different 
properties of these two types of spicules, their differences from observed spicules 
and what needs to be improved in the models. 



1. Numerical Methods and description of the model 

The nature of spicules observed at the solar limb have long been a mystery, 
in this paper we discuss two types of jets that occur naturally in 3D num- 
berical models of the solar atmosphere. The MHD equations are solved in a 
model spanning the upper convection and corona using the Oslo Stagger Code 
(OSC). In addition, this code solves a rather realistic NLTE radiative transfer, 
in cluding scattering, and therm al conduction along the field lines as explained 
in iMartinez-Svkora et all (|2008l ). 



The models described below have a grid size of 256 x 128 x 160 points 
spanning 16 x 8 x 16 Mm 3 . The grid is uniform in the horizontal direction 
with a grid spacing of 65 km. In the vertical direction the grid is non-uniform, 
ensuring that the vertical resolution is good enough to resolve the photosphere 
and transition region with a grid spacing of 32.5 km, while becoming larger at 
coronal heights. At these resolutions the models have been run for roughly 1.5 
hours solar time. We have created two models with different average unsigned 
field in the photosphere; one with 16 G (A2) and the other 160 G (Bl). In 
addition to this ambient field, we in troduce a magnetic flux tube into the lower 
boundary at the bottom boundary ( Martmez-Svkora et al. 20081 ). 



2. Results and discussions 



In the models we found two types of spicule-li ke structures, i.e. the so-called 
type I (jMartmez-Svkora et al1l2009l ) and type n (jMcIntosh et al.1 120071 ). 

A synthetic image of Ca n at the limb is shown in fig [1] which shows the 
two types of spicules. These structures look rather similar to what is observed 
at the limb in the Sun in Ca II. 
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Figure 1. Synthetic image of Ca n H from the limb of the model. Observe 
the two types of spicules, type I located at x = 14 Mm and type II at x = 
7 Mm. The synthetic image is done with MULTL3D. 



Table [T] shows the differences between the two types of the spicules in our 
models and compared to observations. Th e reader is referred to work that ha s 
recently been completed related to spicules; iRouppe van der Voort et a l. (2007); 
De Pontieu et al.1 (|2007l h lHansteen et al.l (|2006l h iMartmez-Svkora et al.l (|2009l ). 



Type I 


Type II 


Observations 


150 examples in both 
models 


2 examples only in Bl 
model 


Type II ubiquitous 


Length w [0.4, 1.5J Mm 


Length ~ 5 Mm 


Type I are longer 


Duration ~ [2,5] min 


Duration ~ 1 min 


Type I have longer dura- 
tions 


Parabolic profile in time 
(deceleration) 


Complex velocity profiles 
due to acceleration at dif- 
ferent height 


Seems to agree (see bibli- 
ography) 


Up-downflow profile 


Only upflow 


Seems to agree 


Velocities ~ [5,35] km/s 


Velocities ~ 150 km/s 


Type I reach larger veloc- 
ities 


Observed in Ca n 


Counterpart in Transi- 
tion region emission lines 


Seems to agree 


Driven by magneto- 
acustic shocks 


Reconnection 


Similar drivers suggested 



Table 1. Properties of the two types of spicules "observed" in the models 
and compared with observations. 



Most likely, the differences between the two types seem to agree with the 
observations. However, a deeper study needs to be done with the type n spicules 
found in the model (work in progress). Moreover, a closer comparison with the 
observations is required. It is interesting to note that the appearance of type I 
does not show a clear preference between models with or without flux emergence, 
while type n only appear in the model with the largest ambient ambient field 
(Bl), and only after emerging flux cross the photosphere. Spicules of both types 
in the models are located at the footpoints of the atmospheric coronal loops, 
where the field lines are open field lines or at least penetrate into the corona. 
Moreover, the footpoint that is closer to the emerging flux tube is the one that 
shows most jets. The type II spicules shows a corresponding nearby hot loop 
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Figure 2. Histograms, normalized to the total number of spicules, for de- 
celerations, maximal velocity, maximum lengths, and duration, from left to 
right, respectively, measured from the two models (Bl dash and A2 dash-dot 
line) and the sum (continuum line). The vertical line is the median value of 
the distribution. The two models show different distributions of the decelera- 
tion, maximum lengths, and duration, as well as some differences in maximum 
velocity. The model Bl shows on average slightly lower decelerations, shorter 
length, and shorter velocities than A2. 



which also seems observed in the Sun (iDe Pontieu et al.1 120091 ) . The hot loop 



(> 10 6 K) can be observed with coronal emission lines. 

In brief, we can summarize the differences between observations and models 
for type I spicules by noting that the upper limits of the d eceleration, length, du- 



ratio n, and maximal velocity are smaller in the models (iMartmez-Sykora et al 



2009). Histograms for deceleration, maximum length, maximal velocity and du- 



ration for the type I from the model are shown in fig [2l These can be compare d 



with the histograms from the observations done by IDe Pontieu et al.l (|2007l ). 
They show an agreement in the lower part to the histograms, and the differ- 
ences betwe en the Bl and A2 s eems similar to the differences of the two regions 



observed by IDe Pontieu et al. ( 20071 ). However, the models do not fit with the 



upper part of the observed histograms. 

In order to improve the models we consider that the resolution of the box 
is important. The chromosphere is poorly resolved numerically and this affects 
the size of the structures of the spicules. In addition, low resolution might 
bring other effects like the diffusion of the shocks (type i) or of the magnetic 
discontinuity (type II). With higher resolution we expect sharper shocks, larger 
range of velocities, better resolved and more frequent spicules. 

In models, it is also important to take time-dependent hydrogen ionization 
into account in the the upper chromosphere. The ionization of hydrogen in the 
solar chromosphere and transition region does not obey LTE, or instantaneous 
statistical equilibrium, as the timescales of ionization and recombination are 
long compared with HD timescales, especially for magneto-acoustic shocks. The 
shock temperatures are higher, and the intershock temperatures are lower, in 
models where time-dependent ionization is considered. This effect wil l likel y 
change the range of parameter of the spicules type I ( Leenaarts et al.ll2007l ). 



Modeling the chromosphere is strongly important to study prop erly the radiativ e 
losses app r oxima tions, NLTE with scattering as has discussed b\ ICarlsson ( 20101 ); 
Leenaartd (|2010T ). 
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Figure 3. Temperature, ion fraction, magnetic field intensity, and ohm, hall 
and ambipolar diffusion calculated as post-processing in a 2D cut of the model, 
from left to right and top to bottom. Observe that Ohm and Hall diffusion 
are rather important in the lower chromosphere and the ambipolar diffusion 
in the upper chromosphere. 

The partial ionization might have other effects on both types of spicules, 
as well. When considering partial ionization we find that ambipolar diffusion, 
Hall diffusion and ohmic diffusion contribute at differing rates throughout the 
chromosphere. The ratio between these three diffusion terms changes from the 
photosphere up to the transition region (see fig [3]). This will possibly have 
important effects in the chromosphere as the parameters controlling reconnection 
and the damping of waves change. 

Finally, we note that the range of ambient magnetic field structures that 
have been modeled only form a small subset of those expected when consid- 
ering supergranulation, plage and the chromospheric network. In addition, a 
continuous weak magnetic flux emergence may need to be added, since it has 
been observed in the models that chromosphere and transition region heights 
are considerably increased with flux emergence. 
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